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A holographic dark energy from Ricci scalar curvature called Ricci dark en- 
ergy was proposed recently. In this model the future event horizon area is 
replaced by the inverse of the Ricci scalar curvature. We study the evolution 
of equation of state of the Ricci dark energy and the transition from decel- 
erated to accelerated expansion of the universe in the Brans-Dicke theory, 
which is a natural extension of general relativity. We find that the current 
acceleration of our universe is well explained. 
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1 Introduction 



Why cosmological constant observed today is so much smaller than the 
Planck scale? This is one of the most important problems in modern physics. 
In history, Einstein first introduced the cosmological constant in his famous 
field equation to achieve a static universe in 1917. After the discovery of 
the Hubble's law, the cosmological constant was no longer needed because 
the universe is expanding. Nowadays, the accelerating cosmic expansion first 
inferred from the observations of distant type la supernovae [T] has strongly 
confirmed by some other independent observations, such as the cosmic mi- 
crowave background radiation (CMBR) [2] and Sloan Digital Sky Survey 
(SDSS) [3J, and the cosmological constant comes back as a simplest candi- 
date to explain the acceleration of the universe in 1990 's. 

Holographic principle [I] regards black holes as the maximally entropic 
objects of a given region and postulates that the maximum entropy inside 
this region behaves non-extensively, growing only as its surface area. Hence 
the number of independent degrees of freedom is bounded by the surface area 
in Planck units, so an effective field theory with UV cutoff A in a box of size 
L is not self consistent, if it does not satisfy the Bekenstein entropy bound j5] 
(LA) 3 < S BH = rrL 2 M^, where M~ 2 = G is the Planck mass and Sbh is the 
entropy of a black hole of radius L which acts as an IR cutoff. Cohen et.al. 
[B] suggested that the total energy in a region of size L should not exceed 
the mass of a black hole of the same size, namely L 3 A 4 < LM 2 . Therefore 

the maximum entropy is S B ^. Under this assumption, Li [7] proposed the 
holographic dark energy as follows 

p A = 3c 2 M p 2 L- 2 (1) 

where c 2 is a dimensionless constant. Since the holographic dark energy with 
Hubble horizon as its IR cutoff does not give an accelerating universe [8], 
Li suggested to use the future event horizon instead of Hubble horizon and 
particle horizon, then this model gives an accelerating universe and is con- 
sistent with current observation [7J [5]. For the recent works on holographic 
dark energy, see ref. [IHl EH E2] • Recently, Gao et.al. [13J suggested that 
the dark energy density may be inversely proportional to the Ricci scalar 
curvature, and they call this model the Ricci dark energy model. They 
have shown that the causality problem disappears in their model and the 
result is phenomenologically viable. Since the Brans-Dicke theory is a nat- 
ural extension of Einstein's general relativity and can pass the experiments 
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from the solar system |14J, it is worthwhile to investigate this model in the 
Brans-Dicke theory. Therefore, in this letter, we study the evolution of the 
equation of state of the Ricci dark energy and the transition from decelerated 
to accelerated expansion of the universe in the Brans-Dicke theory. For the 
relevant works on holographic dark energy in Brans-Dicke theory, see ref. 

nanaini nana EDi i2i]. 

This paper is organized as follows. In section 2 we briefly review the 
Brans-Dicke cosmology. In section 3 we study the evolution of the equation 
of state of the Ricci dark energy and the transition from decelerated to ac- 
celerated expansion of the universe. In the final section we will give some 
discussions and conclusions. 



2 Review on Brans-Dicke cosmology 

The Brans-Dicke theory of gravity is a natural extension of Einstein's general 
gravity. The action for the Brans-Dicke theory with a perfect fluid [22] in 
the Jordan frame is 



S = Id Xy/g 



V <f>V M <f> 
ftR-a; " V )+£ 



16tt V $ 



M 



(2) 



where $ is the Brans-Dicke scalar field representing the inverse of Newton's 
constant which is allowed to vary with space and time and u is the generic 
dimensionless parameter of the Brans-Dicke theory. The Lagrangian Cm rep- 
resents the perfect fluid matter. In the Jordan frame, the matter minimally 
couples to the metric and there is no interaction between the scalar field $ 
and the matter field [15] . The equations of motion for the metric g^ v and the 
Brans-Dicke scalar field $ are 

1 8vr 
2 9 ^ R = "¥ 



f~1 __ D _„ D_ rpM I O rpBD 



Svr ( 3 ) 
V, ; V M $ = T M » 
2cu + 3 " 

Here the stress-energy- momentum tensor for the matter are defined as usual, 
namely = (2/^g)5(^gC M )/Sg^. The explicit form of Tf£ is 

= (Pm + PhiW^Uy + pm9hv (4) 

where pu and pm denotes the energy density and pressure of the matter 
respectively and £L is a four velocity vector normalized as U^lf 1 = —1 . 
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However, the tress-energy-momentum tensor for the Brans-Dicke scalar field 
$ can not be defined in a similar way due to the non-minimal coupling term 
y^$i? in the action which obstructs the separation of the scalar $ Lagrangian 
from the tensor Lagrangian. The explicit form of T^ D is 



rpBD _ J_ 

^ " 8tt 



(5) 

Note that the Einstein' general relativity will be recovered in the u — ► oo 
limit of the Brans-Dicke theory. 

Assuming our universe is homogeneous and isotropic on large scale, we 
work with the Friedmann- Robertson- Walker (FRW) spacetime 



ds 2 = -dt 2 + a 2 (t) 



<h ' 2 +r 2 (dd 2 + stn 2 dd<p 2 ) 



kr 2 



(6) 



where a(t) is the scale factor of the universe. Here k = —1, 0, 1 represent that 
the universe is open, flat, closed respectively. For simplicity we assuming 
$ = $(£), then in the FRW spaceime the field equations take the forms 
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<M UJ I <3> \ 87T k 



H — \ — — \ — \ =— PM -- (7) 



$ / 6 U j 3$ rJ " a? 

® + 3H$ = {p M - 3p M ) (8) 

p M + 3H(p M +p M ) = (9) 

where H = a/a represents the Hubble parameter and the overdot stands for 
the derivative with respect to the cosmic time t. The first equation (J7J) cor- 
responds to the Friedmann equation, the second equation (jHJ) is the equation 
of motion of the Brans-Dicke scalar field. The last is the conservation law 
WT^ = for the matter fluid. We assume the Brans-Dicke scalar field can 
be regarded as a perfect fluid with stress-energy-momentum tensor 

T* D = (pbd+Pbd)U^U^+Pbd9„ = T* U D /G 
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where G denotes the value of Newton's constant and its energy density and 
pressure are given by 



Pbd = 



Pbd = 




Finally, the Bianchi identity V^G^ = , which plays the consistency relation 
[T6l fm [22], leads to an energy transfer between the Brans-Dicke field and 
other matter. 

p BD + 3H(p BD + p BD ) = ~^qPm ■ (10) 
Any physical solutions should satisfy these equations (17|)- (TTU|) simultaneously. 



3 Ricci dark energy in BD theory 

Recently, Gao et.al [33] have proposed a holographic dark energy model in 
which the future event horizon is replaced by the inverse of the Ricci scalar 
curvature. This model does not only avoid the causality problem and is 
phenomenologically viable, but also solve the coincidence problem of dark 
energy. The Ricci curvature of FRW universe is given by 

R = -6(H + 2H 2 + A) . (ii) 

They introduced a holographic dark energy proportional to the Ricci scalar 

where the dimensionless coefficient a will be determined by observations and 
they call this model the Ricci dark energy model. Solving the Friedmann 
equation they find the result 



= ^PA = T^-fi m oe- 3 * + /oe- {4 -^ (13) 
6Hq I — a 
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where Q m Q = 87rGp m0 /3-ffo' x = \na and fo is an integration constant. 
Taking the observation values of parameters they find the a ~ 0.46 and 
f ~ 0.65. The evolution of the equation of state w\ = p\/ p\ of dark energy 
is the following. At high redshifts the value of w\ is closed to zero, namely the 
dark energy behaves like the cold dark matter, and n,owadays wa approaches 
— 1 as required and in the future the dark energy will be phantom. Further 
more this model can avoid the age problem and the causality problem. 

It is worthwhile to investigate the Ricci dark energy in the framework of 
the Brans-Dicke theory. Note that in Brans-Dicke theory, the scalar field $ 
plays the role of Newton's constant ($ ~ 1/G), so it is natural to modify 
Ricci dark energy (TT2"j) in the Brans-Dicke theory as the following 



The energy density of the perfect fluid matter pu in ([7]) contains the cold 
dark matter p m , radiation p 7 and the Ricci dark energy p^. Every component 
itself satisfies the conservation law (J9]). Let us assume the Brans-Dicke field 
vary with time as a power law of the scale factor [20] , and let x = In a 

$(t) = % e nx = ^e nx (15) 

Gr 

Here we set the present scale factor a = 1, so $ — 1/G. One can see that 
in order to be consistent with observations on Newton constant, the value of 
n should be small. Then eq.([7]) becomes 

H\l + n - V) = ^p M e~ nx - ke~^ (16) 
6 6 

and the Ricci dark energy can be rewritten as 

PA = ^ Q(# 2 )' + 2# 2 + ke~^ e» , (17) 

where prime denotes the derivative with respect to x. The energy density of 
cold dark matter and radiation are solved by conservation laws 

p m = pmoe~ 3x , p 7 = p y0 e~ Ax (18) 
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where p m0 and p y0 are the present energy density of matter and radiation 
respectively. Putting eq. (fT7j) . f[T5|) in eq. flTEl) . we get 



H 2 (l +n- |n 2 ) = a Q(# 2 )' + 2H 2 ) + (a - l)ke 



2x 

(19) 

+ Hi (fi m0 e- (3+n) * + n j0 e~^) 
where H is the present Hubble parameter. Solving the eq. ffT9l) we get 

r2 = a-]_ ke _ 2x + 2_Hl n e -(3+n )x 

A - a (-a + 2A + an) 

2-/^/~n ( a \ \ r> 2(2a — A)x 

where A = 1 + n — ^n 2 and go is an integration constant determined later. 
Then the energy density of the Ricci dark energy is 



_ 87rG a(l-w) 3x an ix 

SHq 2A + a(n — 1) 2A + an 



A — a 



(21) 



where Q^o = h/H^. Using the conservation law, 

Pa + 3(pa+Pa) = (22) 
we get the pressure of the dark energy 



8nG an 4x a(n + l){A-l) {n _ 2)x 
llf PA = -2ATa^^° e A^a Qk ° e 



— n 



(23) 



Notice that the equation of motion of the Brans-Dicke scalar field ([8]) can 
be also rewritten in terms of x 

*V + *,) + In^Y = + §g, A - *gp A ) (24) 



Putting eq. (l2"Ui) . (|2"T|) and (I2"3"j) in eq. (l2"lj) . we obtain the equation for n 
i c\ ( OL — l 3a A — 1 \ _ 2r 

+ fl„(„ + 3) - -**-) ™^ +n) \ + „(H + 1 ) ^ (25 ) 

V2 v ; 2cj + 3; 2A + a(n- 1) V2 / 2,4 + an v ; 



( n(n + 1 + -) - ^(n + —)) g^ 3 ** = 
a 2c<j + 3 ay 

Local astronomical experiments set a very high lower bound on u [23J The 
Cassini experiment [H] implies that u > 10 4 . Likewise, a slow fractional 
variation of $ will lead to a small fraction variation of G, consistent with 
observations. Therefore, the interesting case is that when |n| is small whereas 
u is large so that the product n 2 uj results of order unity [19], namely in the 
limit of io — > oo and n 2 oj ~ O(l). So A ~ 1 — \n 2 u and eq. fl25|) becomes 

2{a-l)Q k0 _ 2x 3Q 



- m —Sa- 
ne + 1— -ne 



1 — 7;n 2 u — a (2 — ^n 2 uj — a) 

+ 7^ — rVr ne +So i + Ue (4 «+ 3- > = 

(2 — ^n z uj) \ a ba J 

We see that as long as n is small enough, i.e. smaller than the exponential 
factor, the above equation is always approximately satisfied. But, if the 
exponential factor becomes very large, the above equation will be no longer 
satisfied, then our results will be invalid. In the following we assume that 
this equation has a exact solution at x = x , 



2(a - l)n ko e~ 2xo 3n m0 e- 3x ° 



1 — ^n 2 uj — a (2 — ^n 2 



n^uo — a 



(27) 



(2 — ±n 2 u) \ a ba J 
We will use this equation to constraint our parameters defined later. 
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One can check that the consistent condition eg. (1101) written in terms of 



x is 



QH 2 (n 2 uj -n + n 2 ) + (Hj(n A u - 3n) 

= 6H*ne- nx (Q A + n m0 e~ 3x + fi^e" 4 *) 
Putting eq. (j2"0"]) and f[2"Tl) in eq. f[2"Bl . we obtain the equation for n 



' — n uj + 5n + n 



a — 1 
A — a 



Qna- 



A—l 

'A- a 



o P - 2x 



(28) 



+ 



+ 



2{-2n 2 uj + 8n + An 2 - n 3 u) 
—a + 2A + an 



YlnA 



2A + a(n- 1) 



2(-3n 2 uJ + lln + 4n 2 - n 3 u) \2nA 



2A + an 



2A + an 



(29) 



+ 



3)n 2 u + n 2 + (11 6A)n 



W" 4+ ^ = 



Taking the same limit as before, we get 



n 2 u(l — a) 

b 



■VtkQe 2x 



4n 2 uj 



2-^-a 



-£l m oe 



-3x 



&n 2 uj _ 4x 2 n 2 uj 2 r _ 4+ 2_nW 

^fi 7 oe + - - 3)n 2 ug e { 4+ « *, = 

2 — ^ « 3a 



(30) 



Since n 2 u is of order unity, the above equation can not be always satisfied, 
so we assume this equation is valid near x — 0, which means it is satisfied 
nowadays, i.e. when a = a = 1 



fcO 



4fi 



ITlO 



6Q 



70 



n w 

l-^-a 2-=^ -a 2-^V" 3a 

O o o 



2 - 2 



3 U> = (31) 



In fact, one shall see that the above equation is really satisfied. Notice that 
this equation can not give constraint on the parameters we defined later, it 
is just used for consistent check. 
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Let's come back to the energy density of the dark energy in eq. (l2ip . 
Taking the same limit and keeping the first order of n, we obtain 

2 — ^n z u — a 2 — ±n 2 uj — a 2 — ^n z uj 

n 2 u)a _ 2 . 1 2 . _ ( A_2 + r£u L) 

n z uo + ba — 6 6 

(32) 

Here we can see that when n is exactly zero and uo is finite, the energy density 
comes back to the result in [T3]. Compare with their result, the contribution 
of radiation on the right hand side of the above equation is only the term 
with the first order of n, which means in the small n limit, the radiation 
part is not important. In another word, radiation contribution to the energy 
density is the high order effect compare to the matter part. The matter terms 
on the right hand side of eq. (I3"2~l) consist of two parts, one for zeroth order of 
n and the other the first order. For simplicity we just keep the terms with 
zeroth order of n, namely 

n A = C n m0 e- 3x + D n kQ e- 2x + F e^ +C "> (33) 

where we have define C = - i % , D Q = 2 n2u i a a and F = o (l — \n 2 uj). 

Now we can rewritten the equation of motion fT2"T|) with these three pa- 
rameters as follows 

4C (1 - D )n ko e- 2x ° + 3C (1 + C )n m0 e-^ + 2C iy 4 *° 
2C - (C - 1)A> 

(34) 

F (l + 3C )e(- 3 + g o' 1 ^ _ 
2C - (C - 1)D 

and the consistent equation ( !3~TT) can be also rewritten in terms of Co, Dq 
and Fq. 



2C n 



;i - D )n k0 - 2(1 + c )n m0 - 3a 



70 



+ 2F (1-2C 



2C - (Co - 1)D C 

where we have used 



(35) 



2C — (C — 1)D 2 , n-i\ n /oe\ 

a = — — , n u = 3 (1 - C ) D (36) 

l + o 
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and 

90 - 2 - (1 - C^)D (3?) 

The pressure of the dark energy in eq. (l23l) can be also rewritten in terms of 
these three constants and in the same limit as before 

8nG n 2 ua „ „ ( 2 n 2 u x / 



Hq ^ A n 2 u> + 6a — 6 



IW-* - g - ^ - l) 9o (l - ^) e('^> 



= -D n k0 e- 2x - C - 1 F e(- 3+c o" 1 )^ 

(38) 

Defining the equation of state of the dark energy w\ = Pa/ Pa, we get 

' [AAtoe- 2 * + C 1 F e^ +c o 1 ) x ] (39) 



30/ 



There are three parameters Co, D and F Q to be determined by eq. fl33|) . (1341) 
and (l39lU n terms of the present values of cosmological parameters as follows 



wao = [Dofifco + C ^q] 

<J"ao (40) 

4(1 - D )n k0 e- 2xo + 3(1 + C )a m0 e- 3xo + 20 y0 e- 4a;o 
+ (Cq 1 + 3)F e { - 3+c ° 1)xo = 

and the consistent equation 

Q k0 — 2f2 m0 — 30 7 o — 3^a Oao — 2Oa = (41) 

Note that when fl k0 ~ 0, f2 70 ~ and Qao ~ 3f2 mo , this consistent equation 
requires % ~ —0.9 which is observational accepted, in other words, this 
consistent equation is really satisfied at present day. Then we get 

Oao — D £l k0 



O m o — 3WA0^A0 — -Do^fco 

(42) 



F = —C (3waoQ\o + A)Ofco) 
11 



Here we can use eg. PS)) to eliminate C and F in the last equation of (l4"0"j) . 
Then we get the last equation as equation of D Q . In fact this equation 
contains another parameters xq, and both Dq and xq determines the zero 
point of the function on the left hand side of this equation. So we have 
the freedom to choice D or D Qko, and solve the value of xq to satisfy the 
equation. Generally, the value of Xq cannot be zero, which means the equation 
is not exactly satisfied presently, but one should not forget that actually we 
have a very small factor n multiplied on the left side of the last equation of 
( 140]) . see eq. (!26|) . So that even though this equation is not exactly satisfied 
at x = 0, i.e., at a = ao — 1, it can be approximately satisfied as long as n is 
small enough. When Q k0 = 0, the last equation in (HUj) does not contain the 
.Do parameters any more, and the value of Xq will change, but again small 
enough n will make the equation approximately satisfied at present. Define 
a function of Co, -Do and F as 

N- l {D n k0 ) = 4(1 - D )n k0 + 3(1 + C )tt m0 + 20,0 + (C^ 1 + 3)F (43) 

Then n should be much less than the value of N, namely n << N, to make 
the equation of motion fl26|) approximately satisfied at x — 0. Note that we 
can use eq. p2|) to eliminate C and F in eq.(T43]). 

As an example we will consider the following values of cosmological pa- 
rameters [2J : n A0 = 0.72, tt m0 = 0.27, fi 70 = 10~ 4 , fi fe0 = 10~ 3 , w A0 = -1. 
We draw a 3D graphic to show the evolution of the equation of state w A q 
with parameter D Qko varying in Fig.l. 
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Figure 1. The evolution of the equation of state w\o with parameter Dottko 
varying. The redshift z is defined as a = (1 + z)~ l . 

Here z = corresponds to the present day. We also choice some typical 
values of D Q k0 to plot the evolution of w\ in Fig. 2 



Equation of State W A (D„n k( , = -0.5 ) Equation of State W A (D fl k0 = -0.3 ) 

0.5 0.5 




-1.5 1 -1.5 

-1 0123456 -1 0123456 

Z Z 
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Equation of State W A (D Sl k0 = 0) Equation of State W A (D il k0 = 0.2 ) 




-1 0123456 -1 0123456 

Z Z 

Figure 2. The evolution of the equation of state w\q with different values 
of parameter D Qk . The dashed line corresponds to w\ — — 1. 

One can see that when D^Qko — 0, the result is the same as that in ref.[T3] 
as follows. The equation of state is nearly zero at high redshifts, so the dark 
energy behaves like matter. w\ approaches — 1 at z ~ 0, and in the future w\ 
will be less than —1. When DqQ^o ~ 0.2 the whole curve is almost the same 
as that when -Do^feo = 0. The interesting things is that when -Do^fco ~ —0.3, 
w\ will never less than —1, and when D Qko ~ —0.5, w\ cross —1 twice. 

As we discussed before the power n in eq.f ll5l) should be much less than 
the function iV defined in (J32D to make the equation of motion satisfied at 
present day. We plot the values of iV with respect to D fl k0 in Fig. 3. One 
can see from the figure that as long as n << 0.1 , which is often the case, n 
is much less than N. 
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The upper bound N of n 




-0.6 -0.4 -0.2 0.0 0.2 

Figure 3. The upper bound N(DoQko) of n with different values of 
parameter D^Qko- The dashed line corresponds to N(DoQko = —0.3). 

In Fig. 4 we plot the evolution of the deceleration parameter defined 

a 



Q = 



a 



H 2 



which can be rewritten by 



H H' 
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Deceleration parameter q 




0.0 0.2 0.4 0.6 0.8 1.0 

Z 



Figure 4. The deceleration parameter q with different values of parameter 
D Qko- The horizontally dashed line corresponds to q = 0. The solid, 
dashed, dotdashed and dotted curves represent D Q k0 = —0.5, —0.3, 0, and 
0.2 respectively. 

From Fig.4, one can see that our universe speeds up as z ~ 0.4 ~ 0.8 for 
different values of DqQ^o- This is consistent with the joint analysis of SNe + 
CMB data z T = 0.52 ~ 0.73 HQ [25]. 

4 Discussion and Conclusion 

We have used different values of D Q ko to study the evolution of equation of 
state of the Ricci dark energy and the transition from decelerated to acceler- 
ated expansion of the universe in the framework of Brans-Dicke theory. The 
parameter DqQ^q should be determined by further observations. When D 
or Qko vanished, we come back to the result in ref . [TS] , which means the role 
of the Brans-Dicke scalar field played can be ignored when spatial curvature 
vanished. Once DoQko does not equal to zero, the evolution of the equation 
of state of the dark energy will be changed. Of course, the large absolute 
value of -D fifco are ruled out by current observation. We also find that the 
contribution of radiation is from the next order of the small quantity n which 
comes from the assumption of the power law of the Brans-Dicke scalar fields, 
so that it can be ignored approximately. 
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